Interferometry using spinor Bose-Einstein condensates 
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We study the time-evolution of an optically trapped spinor Bose-Einstein condensate under the 
influence of a dominating magnetic bias field in the z-direction, and a perpendicular smaller field 
that couples the spinor states. We show that if the bias field depends quadratically on time, the 
relative phases of the spinor components affect the populations of the final state. This allows one to 
measure the differences in the time-evolution of the relative phases, thereby realizing a multi-arm 
interferometer in a spinor Bose-Einstein condensate. 
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I. INTRODUCTION 

The realization of Bose-Einstein condensation (BEC) 
in optically trapped ultracold atomic gases [l| has made 
it possible to study the spinor dynamics of multicompo- 
nent, or spinor, BECs. The theoretical framework de- 
scribing their dynamics provides a set of coupled time- 
dependent Gross-Pitaevskii equations for the different 
components of the order parameter @, S 3- The in- 
ternal dynamics can be studied experimentally by taking 
advantage of the coupling between the hyperfine spin of 
the condensed atoms and an external magnetic field. By 
choosing appropriately the magnetic field it is possible 
to control and detect the state of the spinor conden- 
sate. The time scale characterizing changes in the ex- 
ternal magnetic field can be made much shorter than the 
time scale of the spin mixing dynamics ^] . Then the lat- 
ter can be neglected, resulting in a simplified theoretical 
description. 

One possible experimental configuration consists of an 
initially strong magnetic field in the z-direction (bias 
field) and a weaker perpendicular field (coupling field) in 
the xy-plane. The spinor states become Zeeman-shifted 
by the bias field, and adjacent states are coupled reso- 
nantly when the bias field actually crosses a zero-value 
point (i.e., reverses its direction). Then the spinor dy- 
namics can be described by the Landau-Zener level cross- 
ing model (LZ) 0,13, i]- Note that this is not the only 
possible scenario, and towards the end of the paper we 
discuss other possibilities. 

In this paper we consider a level crossing model known 
as the parabolic model, in which the energy levels have 
a quadratic time-dependence. Depending on the value 
of the parameters the energy levels either do not cross, 
cross at one point, or cross twice. In the double crossing 
case the propagator can be, to a very good approxima- 
tion, obtained by applying the LZ model once at each 
crossing and taking into account the dynamical phase 
accumulated by the different spinor components between 
the crossings |9| . The standard version of the parabolic 
model involves two energy levels. However, we are inter- 
ested in applying it in the context of spinor Bose-Einstein 



condensates, which have an odd number of energy levels. 
We thus generalize the parabolic model to these situa- 
tions using the Majorana representation (gI. IToI. [ill. IT3| . 

The new feature we want to point out is that by using 
the parabolic model we can realize a multi-arm interfer- 
ometer in a spinor Bose-Einstein condensate. The role of 
the two crossings is to mix spinor components, thus they 
work as beam splitters in quantum optics p^ . Between 
the crossings the spinor state changes due to the exter- 
nal magnetic field and the time-evolution determined by 
the Gross-Pitaevskii equations. The original parabolic 
model does not take into account interparticle interac- 
tions or other additional phase effects, so in this paper 
we analyze their role. 

The paper is organized as follows. In Section [IT] we 
introduce the Majorana representation and apply it to 
the Landau-Zener model. In Section [ni] we describe 
the two-level parabolic model and its generalization to 
higher dimensions. Finally in Section IIVI we show how 
the parabolic model can be used to detect the phase dif- 
ferences of the spinor components developed between the 
crossings. 



II. THE MAJORANA REPRESENTATION AND 
THE LANDAU-ZENER MODEL 

Any two-level system Hamiltonian, with explicit time 
dependence, can be written as 



(2) 
i I 



(1) 



where 7*^^^ is the identity matrix and S"!^^ = ai/2 is the 
spin- 1/2 operator in the ith direction, given in terms of 
the Pauli matrices Ci. The time evolution determined 
by this Hamiltonian is obtained by solving the time- 
dependent Schrodinger equation 



ih 



dm 

dt 



(2) 



The solution can be written in terms of an evolution oper- 
ator, or propagator U^^''{t, to), which generates the state 
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at time t when the state at a previous time to is given. 
The general unitary form of the propagator for a two- 
level system is 



a P 
-P* a* 



= 1. 



(3) 



Its exact or approximate expression is known only for a 
few specific cases of Hamiltonian ([IJ . When these mod- 
els are generalized into a higher number of levels forming 
the spinor components, the elements of the correspond- 
ing propagator can be given in terms of the elements of 
t/*^^-' using the Majorana representation. Explicitly, if the 
Hamiltonian of an n-level spinor system is 
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An) 



(4) 



with A and Di as in Eq. ([T]), and S'^"'' are the spin S 
(S — 2n + 1) angular momentum operators, then the 
elements of the propagator [/'^")(t,to) are given by 



fc=0 1=1 

\/(n - - l)!(fc + l)\{n - k - 2)\ 



-i~i 



km{n - i - k)l{i - I - 1)1 



(5) 



where we have not explicitly denoted the time- 
dependence. This result has been derived in [1, (see 
also [nl, ) . We emphasize that Eq. ([5]) is valid only if 
the Hamiltonian of an n-level system can be written as 
in Eq. This is true for bosonic atoms in an external 
magnetic field as long as we have only the linear Zeeman 
shift. 

Let us now apply the Majorana representation to the 
Landau-Zener model [1, 0, Hi, i.e., the hnear crossing 
model. The LZ Hamiltonian is 



4'i(i) 



Xt 



Vo 
-Xt 
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where A and Vq are positive constants. Many different 
approaches have been developed to study the dynam- 
ics of such a two-level system. In this paper we work in 
the instantaneous eigenstate basis, i.e. the adiabatic basis 
of Hamiltonian ([6]), since it allows us to use the results 
of Ref. 9] directly. In the LZ model as well as in the 
parabolic double crossing model the state vectors of the 
adiabatic basis coincide (apart from possible sign differ- 
ences) with those of the original diabatic basis of Eq. ^ 
far from the crossing region. If we choose to = — oo and 
t = +00, the propagator in the adiabatic basis is 



f/fi(+oo,-oo) 



where 
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and A is the Landau-Zener parameter. 
Zener phase (p appearing in Eq. ([7|) is 



The Landau- 



TT A, /A\ 



r(i 
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where F denotes the Euler Gamma function. The value 
of (f) decreases monotonically from 7r/4 to as A goes 
from to infinity. Note that the propagator ^ is given 
in the interaction picture (e.g., Eq. (5) in Ref. P). In 
the Landau-Zener model (in both diabatic and adiabatic 
basis) the energy level separations are initially and fi- 
nally infinite, which means that it is convenient to use 
the interaction picture. 

Applying Eq. ^ we can now evaluate the propagator 
for any number n of levels. In the simplest case of n = 3 
we have 



;7fi(+(^,-oo) 



(1 - _R2)e 



ij2 



i(l> 



R^ 
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In the following section we consider the special case of 
the parabolic model in which two consecutive crossings 
appear. We work under the independent crossing ap- 
proximation (ICA), which states that the two crossings 
can be treated as separate linear crossings. The relevant 
time scale in a single linear crossing is given by the Zener 
time tz, 



tz = 



^, A»l, 



h A«l, 



(12) 



where the upper line corresponds to the adiabatic limit 
and the lower one to the sudden limit. Then the condition 
for ICA reads simply as 



tz < tc 



(13) 



where tc is the time interval between the crossings |l6| . 
The parameters of the parabolic model will be chosen to 
satisfy this requirement (see Fig.[T]). 



III. PARABOLIC MODEL 

We now discuss the two-level parabolic model [13, 
[Tsl . [l9| and its extensions to a higher number of levels. 
The Hamiltonian of the parabolic model is 



H 



(2) 

Para 



at^ 



V 

-at^ 



+ b 



(14) 



where a and v are positive and 6 is a real number. The 
sign of b determines whether we have a double crossing 
(& > 0), a single crossing (6 = 0), or a tunneling case 
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(& < 0). In this paper we will consider only the 5 > case 
(Fig.[T]). In terms of the spin-1/2 operators the parabolic 



Hamiltonian can be written as H 

j(2) 



(2) 



2v 5^'^ ■ 



b)Sz . It is convenient to scale the system by defining 
a dimensionless time t ~ vt/h which leads to the new 
Hamiltonian 



H 



(2) 



er — /i 
I 



(15) 



where e — lT?a/v^ and fi = b/v are dimensionless param- 
eters [i]. 

(2) 

Analytic solutions for the propagator UdJ.^„(+oo, —oo) 
have been presented in the literature 0, Il9l | in some lim- 
iting cases. Reference @ provides an analytic solution 
in the adiabatic limit, using the interaction picture and 
the adiabatic basis representation. If generalized, this 
solution suggests that more generally, if condition ([1^)) is 
satisfied, so that the two resonances can be treated as in- 
dependent linear crossings (ICA), their contribution can 
be described by the Landau-Zener propagator of Eq. ([7]). 

— (2) 

Linearizing Hp,^^^ at the crossings we find that the as- 
sociated Zcner time is given in terms of the scaled time 
as 



Tz = 



(16) 



with the upper line corresponding to the adiabatic and 
the lower line to the sudden limit. The distance between 
the crossings is Tc — 2-^/ ji/e. Thus the ICA condition ([13]) 
is satisfied when /i ^ 1 (adiabatic limit) or (n^/eY^^ ^ 1 
(sudden limit). 




FIG. 1: The two-level parabolic model. Solid lines are the 
diabatic energy levels as a function of time, and dashed lines 
are the adiabatic energy levels (featuring an avoided cross- 
ings). We choose the scaled coupling fe > in order to have 
a double crossing case. The Zener time tz and the crossing 
separation tc are also shown. 



We assume that between the crossings only the phase 
evolves. The propagator relative to this evolution reads 



(2) 



1(t/2 





gicr/2 



(17) 



where a characterizes the adiabatic dynamical phase. As 
it is found out in Q and confirmed by our numerical 
calculations, the propagator obtained in this way works 
very well when ICA is satisfied. The total propagator in 
the interaction picture can be written as 



ttT{2)jj{2)jj(2) 
^ LZ ^-ph ^ LZ 



a 

-(3* 



a* 



(18) 

where u]^^^ indicates the transpose of ([7]). The reason 
of using the transpose lies in the fact that after the first 
crossing one of the adiabatic states changes sign. Thus 
the propagator at the second crossing needs to be written 
in the new adiabatic basis. We need to change the signs 
of the non-diagonal elements of (O, or equivalently in 
our case, to take its transpose. 
In Eq. (UHl) we have 



a = e'^^^iR^ + e-''^''+^'l'\l - R% (19) 

P = 2iRy/l-R^ sin((/) + cr/2). (20) 

Here R and (p are obtained from Eqs. ^ and PH]) by 
defining 



A 



1 



2^ 



R = exp 



4^ 



(21) 



Additionally, the dynamical phase is obtained as 




arcsm 



arcsm 



fi ^ — I 



^ — i fi + i 



(22) 



This is the full phase difference accumulated by the 
adiabatic states between the two crossings. Here EE 
stands for EUipticE and EE for EUipticF. These are 

defined as EllipticE[Q;, m] = dd \/\ — m svc? 9 and 

EUipticF [a, m] = dO 1/ \/l - msin^ 9. Equation ^ 
shows that the dynamical phase a can be written as a 
product of Tc and a function which only depends on fi. 
When /i is large we can approximate a « 8//'^/^/ (B^/e) = 
4Tc/i/3, which correspond to the value of a obtained in 
the absence of coupling. Therefore, when the energy sep- 
aration between the two levels at < = is large compared 
with the coupling, it is possible to calculate the dynam- 
ical phase in the diabatic basis. This approximation has 
been used in Ref. [IQi] in the derivation of an approximate 
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propagator for the parabolic model. However, when the 
exact result ([2^ is used, the resulting propagator works 
well also when v is not small compared with the energy 
separation e.t t — 0. 

As we noted with the LZ model, working in the inter- 
action picture is a necessity due to the otherwise infinite 
dynamical phases. If we start either the LZ model or 
the parabolic model in a superposition of states in the 
Schrodinger picture, problems will arise. Since we are 
modelling here the idea of an interferometer, where the 
first beamsplitter creates the superposition, we can avoid 
the infinite phase problem by restricting our study to sit- 
uations where only a single eigenstate on the system is 
populated initially. 

In a two-level model, starting from level 1, the transi- 
tion probability to level 2 is now 



p(2) _ 



We note that the expression for Pi^2 given in Ref. [§| 
contains misprints, the correct form is given by Eq. (j23p . 

As an example we discuss the generalization of the 
parabolic model to three-level systems. The Hamiltonian 
is now 



(23) 



H 



(3) 
Para 



2{at^ - b) V^v 
V2v V2v 
V2v -2{at^- 



(24) 



Using Eq. ([5]) together with Eq. (fT8)) we find the solution 
for a three-level system propagator as 



(3) 



a 

-V2a/3* 
/3*2 



V2a(3 



(25) 



where a and /3 are given by Eqs. and (PO)) . Starting 
again with level 1 initially populated we find that the 
transition probability to level 3 is 



P. 



(3) 



1^3 



16R'^{1~ R^)^smHa/2 + ^). (26) 



If we ignore the rather small A-dependence of (j), we see 
that the parameter A controls the amplitude, and the 
oscillations arise only from a. 



IV. INTERFEROMETRY USING SPINOR 
CONDENSATES 

A. Theoretical analysis 

The order-parameter of a spin-i^ condensate can be 
written as 

F 

E Vw(r,t)|i^,m), (27) 

mF — — F 



with the normalization X]mj=.=-F = "-(r,*), 

where n(r, t) is the total particle density. We assume 
that the condensate is confined in an optical dipole trap 
so that all the components of the hyperfine spin can be 
trapped simultaneously, and they are degenerate in the 
absence of any magnetic fields. In the following we con- 
centrate on an = 1 condensate. In the degenerate 
case the time-evolution of an = 1 condensate is given 
by the following set of time-dependent Gross-Pitaevskii 
equations [lIllHl 



ih 

ih 
ih 



dt 

dt 
dt 



CiPl + \a (V'oV'll + IV'llVl + IV-olVl 
-IV'-llVl), 

C^O + Aa (2^lV'-l^o* + 1^-1 1 Vo + iV'il'^o) , 

C^'-l + Aa {i>lri + IV'-llV-l + I'/'OI'V'-I 

-m^^-i), (28) 



with C = -ft2v2/(2m) + C/(r,<) + Xsn{r,t), = 
4:Trlf{ao + 2a2)/(3m) and Aa = ATTh^{a2 - ao)/(3m). 
The afe's are the s-wave scattering lengths in the scat- 
tering channels with total spin k. The trapping potential 
is denoted by U , and it affects all spinor components 
equivalently. The first term inside each parentheses de- 
scribes spin mixing processes where atoms in to^^ = -fl 
and mp — ~1 states collide and create two atoms in 
mp = state, or vice versa. The number of particles 
in the ith component is by definition Ni = J cfir\ipi\'^. 
From Eqs. we get 



dNi _ dN^i _ 
"df ~ dt ~ 

It follows that 

dNo 



IdNo 
'2~dr 



2Xa 

h 



d-^rIm[V'oV'-iV't 



dt 



< 7^0, 7 = r- 



(29) 
(30) 



where rimax is the maximum density. Using the scatter- 
ing lengths relevant for ^^Rb 21] and a maximum density 
nmax = 10^^ cm~^ we find 7 ~ 90 s^^. If for example 
we have At = 100 /zs the change in the populations is 
less than 1 %. The time-evolution given by the Gross- 
Pitaevskii equations (pS)) is then restricted to changes in 
the phases of the spinor components. Thus we can de- 
scribe it defining the following effective propagator 



U. 



(3) 
GP 



1 









(31) 



where we set one of the three independent phases to 
zero. The phases 9i and 6-1 can be calculated by solving 
Eqs. (EHl). 
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B. The spinor interferometer 

We base our interferometer proposal for the analogy 
between optical beamsplitters and the level crossings. 
The initial single-state signal is split into two parts at the 
splitter-crossing. The two parts travel along the different 
arms of the interferometer, and combine again in the sec- 
ond splitter-crossing, to form a detectable signal. Thus 
any difference in phase evolution along the two paths is 
mapped into variation of populations. The closest coun- 
terpart of our spinor interferometer is the quantum opti- 
cal Mach-Zehnder interferometer. In our case, however, 
the interferometer can have more than two arms. 

Let us consider a few points in our analogue. The 
splitting at the level crossings can be controlled by the 
speed of change for the bias field and by the coupling field 
strength. Thus we are not limited to the 50-50 splitting. 
Here the controlling parameter is the Landau-Zcner term 
R. As for the phase evolution, Eq. (|26|) demonstrates that 
it decouples from the splitting, although we must remem- 
ber that the standard phase evolution a will depend on 
the same parameters as R. On the other hand, we have 
altogether three parameters [a, b and v in Eq. p^ ]. so it 
is possible to fix R and a independently. 

The time-evolution of spinor states between the cross- 
ings represents the internal part of the spinor interferom- 
eter. This is the time region during which the phenomena 
to be investigated take place. As discussed earlier, we 
consider temporal crossing separations which are much 
shorter than the typical time scales of spin mixing pro- 
cesses. Thus we suggest that the contribution ([51]) is 
included to the total propagator, which becomes now 



a 



(3) 

TOT 



(32) 



("3") 

where U^fJ is obtained from ([TT]) through Eq. ([5]). Thus 
the solution is obtained as a sequence of separate events. 

It is important to notice that the propagator p8)) is 
written in the adiabatic basis while (|3ip is written in the 
diabatic basis [following Eqs. (|^]. But the result ([5^ is 
still well-defined. The reason lies in the ICA approxima- 
tion. Since v <C 6, the adiabatic and diabatic state evolu- 
tion coincide approximately outside the crossing regions. 
The only problem arises in the vicinity of the crossings, 
but these time intervals are much shorter {tz ^ tc) and 
do not contribute much to the phase evolution. 

As an example let us start with the state |1, 1) initially 
populated. The final population of the state |1, — 1) is 

-PH-1 = 16i?*(l - RY {sin*x-hcos2xsin2^'} , (33) 
where 



^-2 



(34) 



The result psp depends on the parameters cr and 0, 
which one can control externally, and on the GP phases 
6i and Thus the measure of the final population of 



the |1, — 1) state (the output of the interferometer) gives 
information on the Gross-Pitajevskii dynamics. We can 
also see that the visibility of the interferometric oscilla- 
tions is maximized when R — 1/ ^/2. 

A simple interferometric measurement would be thus 
to set R = l/\/2 to maximize the signal, then to vary x 
(keeping R constant). As one can see from Eq. the 
difference between minima and maxima gives informa- 
tion about 5*. Interestingly, possible fluctuations in the 
magnetic fields affect the Zeeman shifts between adjacent 
spinor states equally, thus affecting only a. 

One should note that because of the sin'^-function, the 
fringes are "sharper" than in the standard two-arm in- 
terferometer. If we go to higher values of spin, such as 
F = 2, the fringes become even sharper. In analogy 
to diffraction patterns, the double slit is replaced by a 
2F + 1-slit. Note that in this case the number of inde- 
pendent relative phases is 2F. 



C. Experimental realization 

In this section we estimate the values of the parameter 
of the interferometer looking at a real situation. The 
effect of the interaction with an external magnetic field 
is described by the general linear Zeeman model 



H^B^ - MbB • F, 



(35) 



where /is is the Bohr magneton, gp is the Lande g-factor 
(equal to 1/2 for ^^Rb and ^^Na), B = {B^,By,B^) is the 
magnetic field, and F = {F^, Fy,Fz) is the spin operator 
of a spin-1 particle (with h = 1). We can identify iJ^*^ 

with H^pl.ra defining — —Av/^Bi By = 0, and B^ ~ 
—A{at^ — b)/^B- With these definitions we get {gp = 1/2) 

., \B.{t^O)\ 



tn = 




(36) 



\Bi 



where \B^\ is the rate of change of Bz{t) at the crossings. 
The Zener time becomes 




B^ 








/ 2h 









(37) 



We estimate the values of e and fi by choosing B^ — 60 
mG and B^ = b ■ 10"* G/s, and B^{Q) = 300 mG. These 
give e « 2 and = 5. The distance between the crossings 
is tc = 24 /is while the Zener time is « 2 /is, which 
shows that the ICA is well satisfied, and our approach 
should be applicable. From the previous values we also 
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get i? « 0.78 which is very close to the value l/\/2 which 
maximizes the interferometric signal in (|33p . 

In our experimental scenario we have considered a con- 
stant coupling field. Alternatively, we could use an rf- 
field scenario in analogy to evaporative cooling [22] . By 
performing the rotating wave approximation and moving 
to a rotaing coordinate frame we get a situation that is 
euivalent to the constant coupling field scenario. Then 
the crossings take place when the bias-field-induced Zee- 
man shifts become resonant with the rf-frequency, and 
the bias field does not have to pass a zero-point value. 
Yet another variation would involve chirping of the rf- 
frequency, as it can also be described as time-dependent 
variation of the energy of the spinor states. 

V. CONCLUSIONS 

In this paper we have presented the basic idea of a 
spinor interferometer. It is based on applying time- 
dependent magnetic fields to provide level crossings 
of spinor states in analogy to optical beam splitters. 
Interferometer-like studies that apply the Landau-Zener 
model twice and show path-related interference are not 
uncommon in physics, see e.g. Ref. [1^ . The spinor inter- 
ferometer is a partial extension into the case of a multi- 



arm interferometer, which is perhaps harder to realize in 
quantum optics. The extension is only partial due to the 
nature of the spinor. For spin-F we get 2F + 1 paths 
and 2F independent phases, but the spinor dynamics is 
nevertheless characterized by the same small number of 
parameters as the corresponding two-state model, as the 
Majorana representation shows. Instead of being a limi- 
tation, one can see this spinor property as an advantage 
when considering possible error sources. 

It should also be clear that one is not limited to the 
parabolic model. We have chosen to use it because then 
we can link the intuitive concept of using the LZ model 
twice within ICA to the full solution of the model in the 
adiabatic limit, as done in Ref. Our next step is to 
consider more carefully how the dynamics of Eqs. ([28]) 
appear in the phase evolution given by Eq. (PT|) . 

Acknowledgments 

The authors acknowledge the Fondazione A. Delia Ric- 
cia and Finnish CIMO (R.V.), EC projects CAMEL and 
EMALI (H.M), and the Academy of Finland projects 
108699 and 115682 (H.M and K.-A.S) for financial sup- 
port. The authors thank Boyan Torosov and Nikolay 
Vitanov for enlightening discussions. 



[1] D. M. Stamper-Kurn, M. R. Andrews, A. P. Chikkatur, 
S. Inouye, H.-J. Miesner, J. Stenger, and W. Ketterle, 
Phys. Rev. Lett. 80, 2027 (1998). 

[2] T. Ohmi and K. Machida, J. Phys. Soc. Jap. 67, 1822 
(1998). 

[3] T.-L. Ho, Phys. Rev. Lett. 81, 742 (1998). 

[4] C. K. Law, H. Pu, and N. P. Bigelow, Phys. Rev. Lett. 

81, 5257 (1998). 
[5] H. Pu, C. K. Law, S. Raghavan, J. H. Eberly, and N. P. 

Bigelow, Phys. Rev. A 60, 1463 (1999). 
[6] E. Majorana, Nuovo Cimento 9, 43 (1932). 
[7] L. D. Landau, Phys. Z. Sowjet Union 2, 46 (1932). 
[8] C. Zener, Proc. R. Soc. Loud. A 137, 696 (1932). 
[9] K.-A. Suominen, Opt. Comm. 93, 126 (1992). 
[10] F. Bloch and I. I. Rabi, Rev. Mod. Phys. 17, 237 (1945). 
[11] N. V. Vitanov and K.-A. Suominen, Phys. Rev. A 56, 

R4377 (1997). 

[12] L. Xia, X. Xu, R. Guo, F. Yang, W. Xiong, J. Li, Q. Ma, 
X. Zhou, H. Guo, and X. Chen, Phys. Rev. A 77, 043622 
(2008). 

[13] C. Gerry and P. Knight, Introductory Quantum Optics 

(Cambridge University Press, Cambridge, 2005). 
[14] A. Meckler, Phys. Rev. Ill, 1447 (1958). 



[15] A. P. Kazantsev, G. I. Surdutovich, and V. P. Yakovlev, 

Mechanical Action of Light on Atoms (World Scientific, 

Singapore, 1990). 
[16] K. Mullen, E. Ben- Jacob, Y. Gefen, and Z. Schuss, Phys. 

Rev. Lett. 62, 2543 (1989). 
[17] V. K. Bykhovskii, F. F. Nikitin, and M. Ya. Ovchin- 

nikova, Sov. Phys. JETP 20, 500 (1965). 
[18] D. S. F. Crothers and J. G. Hughes, J. Phys. B: Atom. 

Molec. Phys. 10, L557 (1977). 
[19] E. Shimshoni and Y. Gefen, Ann. Phys. (N.Y) 210, 16 

(1991). 

[20] T. Isoshima, M. Nakahara, T. Ohmi, and K. Machida, 

Phys. Rev. A 61, 063610 (2000). 
[21] E. G. M. van Kempen, S. J. J. M. F. Kokkelmans, D. J. 

Heinzen, and B. J. Verhaar, Phys. Rev. Lett. 88, 093201 

(2002). 

[22] C. J. Pethick and H. Smith, Bose-Einstein Condensation 
in Dilute Gases, 2nd ed. (Cambridge University Press, 
Cambridge, 2008). 

[23] M. Sillanpaa, T. Lehtinen, A. Paila, Y. Makhlin, and P. 
Hakonen, Phys. Rev. Lett. 96, 187002 (2006). 



